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ABSTRACT 

Employing a catalog of 175 extrasolar planets (exoplanets) detected by the Doppler- 
shift method, we constructed the independent and coupled mass-period functions. It is 
the first time in this field that the selection effect is considered in the coupled mass- 
period functions. Our results are consistent with those in Tabachnik and Tremaine 
(2002) with the major differences that we obtain a flatter mass function but a steeper 
period function. Moreover, our coupled mass-period functions show that about 2.5 
percent of stars would have a planet with mass between Earth Mass and Neptune Mass, 
and about 3 percent of stars would have a planet with mass between Neptune Mass and 
Jupiter Mass (Mj). 

Key words: planetary systems, methods: data analysis, methods: numerical, methods: sta- 
tistical 



1. Introduction 

The mass (size) function, i.e. the differential form for the number of objects as a function of 
mass (size), 

f ( \ dN m 
f M {m) = — , (1) 

am 

has been an important physical property to be investigated in many fields of astronomy such as 
galaxies, stars, asteroids and also dust grains. The importance lies in the meaning that this function 
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is related with the formation and evolution of particular types of objects. Due to the fact that the 
mass function can be studied either through observational techniques or theoretical calculations, 
numerous research projects have been done on this subject. 

For stars, the initial mass function (IMF) is the distribution of stellar masses from one star 

formation event in a given volume of space. Although the star-forming conditions vary with the 

environment, the measured IMF appears to be universal and can be modeled by a power-law, 

dN 

f M (m) = — = c,m- a % (2) 

where c* is a normalization constant, and a* = 2.35 for the well-known Salpeter IMF. According 
to Kroupa, Tout, and Gilmore (1993), 

fuijn) = cim -4 - 5 for m > 1M , 

f M (m) = c 2 m- 2 - 2 for O.5M < m < 1M , (3) 

/ M (m) = c 3 m -L2 for m < O.5M , 
where m is the star's mass, and ci, c 2 , C3 are constants to be determined by the total number of 
stars in the considered system. 

Moreover, there is a new development in astronomy that more than 300 planets are detected 
around solar type stars. The discovery has led to a new era of the study of planetary systems 
and thus triggered many interesting or controversial results of theoretical works (Jiang & Ip 2001, 
Kinoshita & Nakai 2001, Armitage et al. 2002, Ji et al. 2003, Boss 2005, Jiang & Yeh 2007, Rice 
et al. 2008, Ji et al. 2009). For example, many discovered exoplanets have extremely short orbital 
periods. It is likely that they are formed at larger radial distances and migrate to the current 
locations later. However, because the migration timescale is too short, the rapid inward Type I 
migration caused by disc-core interactions poses a serious issue. 

At the time when there were only about 70 detected exoplanets, Tabachnik & Tremaine (2002) 
first used the maximum likelihood method to determine the mass and period functions with the 
assumption of two independent power laws. This work is, in fact, the only one that takes into account 
the selection effect and has intentions to obtain the fundamental mass-period functions. We have 
to note that different from the definition of stellar IMF, the planetary functions are constructed 
through the data of exoplanets from all different systems. 

Without considering the selection effect, Zucker & Mazeh (2002) calculated the linear corre- 
lation coefficient between mass and period and concluded that the mass-period correlation exists. 



- 3- 



The study of correlations within sub-groups in the cluster analysis done by Jiang et al. (2006) and 
Marchi (2007) also confirms the mass-period correlation. Thus, strictly speaking, the mass and 
period functions shall not be two independent power laws. Motivated by the above results, Jiang et 
al. (2007) employed an algorithm to generate a pair of positively correlated /9-distributed random 
variables in order to construct a coupled mass-period function, and was the first one in the field of 
exoplanets to include the correlation into the construction of mass-period functions. 

Moreover, through a nonparametric approach, Jiang et al. (2009) further constructed new 
coupled mass-period functions from 279 exoplanets and presented two main statistic results: (1) 
confirming the deficit of massive close-in planets; (2) discovering that more massive planets have 
larger ranges of possible semimajor axes. However, due to the fact that the selection effect is not 
considered in the above study, it is unclear how strong the statement (2) is. Yeh et al. (2009) 
argued that the planets larger than 1 Mj should be all within the detection limit, and thus implied 
that (2) should be a statistically valid statement. 

Therefore, in order to address the above problem, it would be a great success if the coupled 
mass-period function can be obtained while the selection effect is considered simultaneously. It 
is this goal that motivates this work. The coupled mass-period function will be constructed by 
the same statistic method as in Jiang et al. (2009) and the selection effect would be considered 
through the similar procedure used in Tabachnik & Tremaine (2002). However, as this work is an 
extension from Tabachnik & Tremaine (2002), power laws would be employed as standard forms of 
mass-period functions. In order to consider the selection effect, only the exoplanets detected by the 
Doppler-shift method will be included as our samples. 

We first constructed a reference-based catalog in §2, and discussed the normalization problem 
in §3. Then, the independent mass-period function was constructed in §4. In §5, mass-period 
correlations was studied. The coupled mass-period function was established in §6, and the fraction 
of stars with planets was discussed in §7. We concluded the paper in §8. 

2. A Reference-Based Catalog 

The most well-known catalog of exoplanets is the one maintained by Jean Schneider (Schneider 
Catalog hereafter), i.e. http://exoplanet.eu/catalog.php. This catalog is updated frequently when 
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new detections are reported in refereed papers or conference proceedings. It is also notable that 
Butler et al. (2006a) updated the orbital solutions and compiled a list of 172 exoplanets. In order to 
extend the work of Tabachnik and Tremaine (2002) with more updated data, we here constructed a 
new exoplanet catalog, in which all samples were discovered by the Doppler-shift method, as listed 
in Appendix A. A major principle to construct this catalog is that all exoplanets in our catalog shall 
be reported as new detections in papers of refereed journals. In order to establish such a reference- 
based catalog, we searched and reviewed many published papers. We intended to make sure that 
these papers did report new discoveries, and to check which observational survey the results belong 
to. Thus, all the references listed in our catalog are the papers which reported new detections. 

In our catalog, the first column is the data-set identity and the 2nd column is the name of 
the observational survey. The 3rd and 4th columns are the reference papers and the papers' corre- 
sponding identities. The 5th column lists the number of exoplanets discovered in that corresponding 
paper. Finally, column 6, 7, and 8 give the name, the projected mass M, and the orbital period P 
of exoplanets. Most planets' M and P are obtained from Schneider Catalog. The values in Butler 
et al. (2006a) are used when they are missing in Schneider Catalog. We have to note that the 
exoplanet HD 154345b in Ref. (E-l) is removed here due to its extraordinarily larger period, 10900 
days. 

Besides, we found that some exoplanets were detected by more than one group around the 
same time, they could be reported as new discoveries in two different papers. We checked these and 
made a list in Table 2. However, they are still repeatedly listed in Appendix A, as our catalog was 
constructed based on published papers. Moreover, some papers studied more than one exoplanets, 
among which only one planet is a new discovery. In those papers, additional planets were included 
for a comparison purpose. We also carefully examined these kinds of papers, so that the reference 
papers in the catalog are exactly the papers which discover those listed exoplanets. 

3. The Normalization Problem 

In addition to the selection effect, the number of observed stars in a survey, N s , is one of the 
key parameters to understand the probability that a star could host an exoplanet with particular 
orbital properties. Unfortunately, we found that the exact numbers of stars in surveys were not 
clearly mentioned normally. In most cases, an approximate number of target stars in a continuous 
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long-term survey might be given, but the exact number of observed stars at that time when a new 
exoplanet detection was detected was usually not stated. For example, in Butler et al. (2003), it 
is written as "The Keck survey includes about 650 main-sequence and subgiant stars and later 
in the same section, it is stated as "200 stars have one or more Keck observations but have been 
subsequently dropped from the program In this case, neither 650 nor 450 can be used as the 
N s here. 

Fortunately, Lineweaver and Grether (2003) obtained the numbers of target stars in surveys 
carefully and listed the results in their Table 4. In fact, they also estimated the number of repeated 
stars in different surveys, and finally determined the total number of target stars to be 1812 in these 
surveys. The numbers of detected planets in the corresponding surveys are also shown in their Table 
3, and the total number is 122. Therefore, we decide to set iV ratio = N/N s = 122/1812 = 0.0673 in 
our paper. This value of N/N s would be used when we determine the value of the normalization 
constant c through the related equations, and also as a way to obtain N s in a particular survey for 
a given number of detected exoplanet N. 

4. Independent Mass-Period Functions 

In this section, the method and results based on the assumption of independent mass and 
period power-law functions will be described. 

4.1. The Method 

The procedures in this section, i.e. the analytical approach, the choice of parameters etc., 
follow exactly that of Tabachnik and Tremaine (2002). Here we describe the method in a self- 
consistent way. However, please note that we simply consider the projected (minimum) mass M, 
which satisfies M = M real sini, where M real is the real physical mass of the exoplanet and i is the 
orbital inclination angle, in all calculations in this paper. That is, the mass M in this paper means 
the minimum mass. The probability, dp, that a single star has an exoplanet with mass M and 
orbital period P in the range, [M, M + dM], [P, P + dP], is given by the product of independent 
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power laws on M and P as: 



dp 



' M ' 


—a 


' P~ 


[M _ 




Po 



dM dP 

~W~p' 



(4) 



where c, a and f3 are constants to be determined, and M = 1.5Mj, P = 90 days. We assume that 
there are N exoplanets in the data set, and let 

Xi = \n( 

y^ln^gj for 1 < i < N 

where Mj, Pi are the mass and orbital period of one particular exoplanet. According to Eq.(2) and 
(6) of Tabachnik and Tremaine (2002), the value of Xi — yi/3 shall satisfy 



Xi - Vi/3 > In 



K 



D 



28.4ms 



-l 



-In 



Mo 
Mi 



1, 

+ -In 

3 



lyr 



(6) 



where K D is the detection limit of the considered survey. For convenience, we define the right hand 
side of the above equation to be 



v D = In 



K 



D 



28.4ms- 



ln 



Mo 
Mi 



lyr. 



(7) 



The value of vp can be determined if the detection limit K D is exactly known. However, a range 
of possible values of Kb is usually stated in a paper because it is related with the condition of 
instruments and weather. Similar to Tabachnik and Tremaine (2002), we set the smallest value of 
Xi — Ui/S, i.e. mini<j<Ar {xi — yi/3}, in the considered data to be the value of vd and obtain the 
most likely K D reflected from the data through the above equation. Further, any value of yi shall 
satisfy 



Vi < In 



(8) 



where P max is the upper limit of orbital periods of detectable systems, which is proportional to the 
duration of the survey, as at least two orbits are required for a reliable detection. For convenience, 
we define 



Ud = In 



(9) 



Similarly, it is not easy to obtain a consistent P max from the published papers. We set the largest 
value of yi in the considered samples to be the value of up and get the most likely P max reflected 
from the data. 

On the other hand, there are further constraints on X{ and y{. Let M max be the maximum value 
of the mass in the data set for a considered survey, and define 



,x T] 



In 



" M ma 
Mo 



(10) 
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so we have a constraint on Xi as 

Xi ^ x max . (H) 

Similarly, let P min be the minimum value of the orbital period in the data set for a considered 
survey, and define 



2/min = ^ 

we have another constraint on yi as 



p . 

1 mm 



(12) 



Vi > Vmin- (13) 

Therefore, in the x — y space, the exoplanet probability dp can be expressed as: 

dp = c[e x }- a [e y }- (3 dxdy, (14) 

and the expected number of exoplanets in the area dxdy in a survey of N s stars can be written as 

n(x, y)dxdy = N s dp = cN s e~ ax e~^ y dxdy. (15) 

Thus, as in Tabachnik and Tremaine (2002), the likelihood function L is expressed as 

L^ in{X ,M-[n( xM] , (16) 

Jd 

where the domain Disv + y/3<x< x max , y niin < y < u with 

v = mhii^Ar {xi - f} 

u = max^^Ar {yi} (17) 
u = mm{u,3(x ms , x - v)} 

Finally, through the maximum likelihood method, the values of c, a, and f3 can be determined by 
the following equations: 

<91nL <91nL B\nL , . 

_ = 0, -^ = 0, and — = 0. (18) 

The first one, d In L/dc = 0, can be solved analytically and yielded an analytical expression of c 
as shown in Eq. (16) of Tabachnik and Tremaine (2002). Using this analytical expression of c, 
d In L/da = and d\nL/d/3 = can be solved numerically. 
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4.2. Individual Surveys 

First, the mass-period function is estimated through the data from one given particular survey. 
The previously mentioned equations are used to obtain the estimates of a, (3, and c. Their error 
bars are calculated as the standard deviations via the bootstrap algorithm with the replication size 
B = 2000 (Jiang et al. 2007). The results are shown in Table 1, and each row is for one particular 
survey. 



Table 1. 



Survey 


Data Set 


N 


a 


P 


c x 10 3 


Mms- 1 ) 


Pmax(yr) 


Lick 


(A) 


7 


-0.401±1.051 


-0.144±0.388 


4.466 ± 4.915 


19.794 


14.296 


Coralie 


(B) 


38 


0.332±0.186 


-0.393±0.101 


2.969 ± 0.372 


21.032 


5.759 


Elodie 


(C) 


14 


-0.473±0.827 


-0.343±0.291 


2.327 ± 5.876 


21.400 


7.921 


HARPS 


(D) 


23 


0.081±0.146 


0.097±0.146 


1.733 ± 4.189 


1.128 


2.844 


N2K 


(E) 


14 


-0.071±0.425 


0.077±0.196 


3.320 ± 4.057 


12.998 


3.507 


Keck 


(F) 


48 


0.131±0.135 


-0.326±0.095 


1.362 ± 0.243 


2.9304 


9.184 


AAPS 


(G) 


23 


0.296±0.262 


-0.558±0.202 


1.778 ± 0.715 


10.318 


8.181 


Other 


(H) 


17 


0.168±0.350 


-0.341±0.192 


3.896 ± 3.154 


28.189 


5.957 


Single 




175 


-0.143±0.040 


-0.124±0.044 


1.316 ± 0.080 


1.128 


14.296 



The data sets and total numbers of exoplanets used in the calculation of maximum likelihood 
method are shown in columns 2 and 3. The resulting values of a and (3, i.e. the exponents of the 
power-law mass-period functions, are in columns 4 and 5. To make it clear, their values with error 
bars are plotted in Fig. 1. Table 1 shows that the Lick results have the largest error bars, due to 
the small numbers of planets. Note that Lick Survey has the smallest number of planets here. The 
reason is that we put the planets in Fischer et al.(2001) into Other Survey because both data came 
from Keck and Lick telescopes are used in that paper. We found that, among these surveys, the 
Coralie and Keck results have smaller error bars as Coralie Survey has 38 and Keck Survey has 48 
planets. HARPS Survey has 23 planets, but the error bars are about as small as those in Coralie 
Survey due to the high precision on the measurement of radial velocities. In general, the absolute 
values of a are less than 0.5, and the values of (3 range from -0.6 to 0.1. They both can be positive 
or negative. 
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The final row is the result when all the above exoplanets are assumed to be discovered by a 
Single Imaginary Survey. However, as shown in Table 2, some planets are included in more than one 
data set. Although the summation of numbers of planets from all data sets is 184, there are actually 
175 planets in total. The a and (3 of the result of Single Imaginary Survey are both negative, and 
the corresponding error bars are very small due to a much larger number of samples. 



Table 2. Repeated Exoplanets 





Planet 


Reference ID 


(1) 


HD102117 


(G-12) and (D-5) 


(2) 


HD 196050 


(B-l) and (G-2) 


(3) 


HD 216437 


(B-l) and (G-2) 


(4) 


HD 52265 


(B-4) and (F-16) 


(5) 


HD 192263 


(B-7) and (F-7) 


(6) 


HD 168443 c 


(B-ll) and (F-10) 


(7) 


HD 33636 


(C-2) and (F-9) 


(8) 


HD 37124 c 


(F-l) and (F-12) 


(9) 


HD 92788 


(B-l) and (H-l) 



4.3. Multiple Surveys 

The analysis of individual surveys is generalized to Multiple Surveys here. In this analysis, 
different domains are considered for different surveys in Eq. (16). Thus, from Eq.(16), the likelihood 
function Lj for the survey j is 

Lj = U^i 1 n j (x itj ,y itj )ex.-p[- / n j (x,y)dxdy\, (19) 

where Nj is the number of discovered planets for the survey j, and 



cN- 

nj (x,y) = cN s e~ ax -^ = j^-e^ 

1 v r; 



'ratio 



Thus, the log-likelihood is given by 



lnL = ^lnL j , (20) 
j'=i 
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where J is total number of considered surveys, and we have J = 8 here. When InL is maximized, 
the estimates for the parameters are obtained as: 

a = -0.099 ±0.055 (21) 
P= -0.129 ±0.040 (22) 
c= (2.423 ±0.154) x KT 3 , (23) 

where the error bars are also estimated via the bootstrap algorithm (Jiang et al. 2007). This set of 
a and (3 with error bars are also shown in Fig. 1. 

In order to show the distributions implied by these results and to be compared with the samples' 
histograms, we take integrations for the function n(x, y), and thus the number of planets with masses 
between M 1 and M 2 and periods between P 1 and P 2 is given by 

/•ln(P 2 /Po) /-MMa/Mo) n r 

N[ Mi ,m 2 }[pi,P2] = / n(x,y)dxdy = / / cNge'^'^dxdy, 

J\n(P 1 /P ) Jln(M 1 /M ) J J 



"ln(P 2 /P ) 




af3 \\Mq) \M J J 1 \P ; 1 " ; ' 

After substituting the result of Single Imaginary Survey into the above equation, we obtain the 
expected number of planets in given mass or period intervals as shown in Fig. 2. With mass 
interval 0.5 Jupiter-mass and period interval 150 days, the crosses connected with dotted lines in 
Fig. 2 show the planetary distributions in mass space (see Fig. 2(a)), and also in period space 
(see Fig. 2(b)). For the comparison, the mass and period histograms of 175 exoplanets and the 
corresponding results derived from the Multiple Surveys are also plotted in Figs. 2(a)-(b). 



5. Mass-Period Correlation Coefficients 

In addition to the mass and period histograms and distributions shown in Fig. 2, the exoplan- 
ets' locations on the mass-period space, i.e. M — P space, are often presented, as the mass and 
orbital period are the most important physical parameters of exoplanets. Fig. 3 is the distribution 
of 175 exoplanets in logarithmic space, x — y space, where x, y are as defined previously. The region 
enclosed by four solid lines is the Domain D, where all 175 planets are included. In order to inves- 
tigate the strength of mass-period correlations, the Spearman rank-order correlation coefficients ps 
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are calculated, as shown in Table 3. 



Table 3. 



SI 11 T VP V 


Data Set 


N 


ys 


Lick 


(A) 


7 


0.143 


Coralie 


(B) 


38 


0.351 


Elodie 


(C) 


14 


0.431 


HARPS 


(D) 


23 


0.564 


N2K 


(E) 


14 


0.552 


Keck 


(F) 


48 


0.532 


AAPS 


(G) 


23 


0.553 


Other 


(H) 


17 


0.542 


Single 




175 


0.540 



In general, the values of ps are larger than 0.5 in five individual surveys and the Single Imaginary 
Survey of overall 175 planets. These results imply a strong positive correlation (Cohen 1988). 
Because there is a mass-period correlation, we shall consider the mass-period coupling and construct 
coupled mass-period functions in this paper. 

On the other hand, as pointed in Jiang et al.(2006) and in Marchi (2007), the correlation of 
each group of exoplanets identified in the clustering analysis might be linked with the physical 
mechanisms more easily. The overall correlations could be more difficult to explain. A recent result 
by Marchi et al. (2009) is a good example that two dominated groups of close-in exoplanets are 
found and can be explained by two physical mechanisms successfully. This is an interesting and 
important topic that we would like to do further investigations in the future. However, the main 
point in this section is to re-confirm that there is indeed a mass-period correlation for our selected 
samples, and thus it is necessary to construct the coupled mass-period functions in this paper. 
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6. Coupled Mass-Period Functions 

We consider the mass-period coupling and construct coupled mass-period functions here. We 
define variables x = ln(M/M ), y = \n(P/P ), and their probability density functions as: 

f X ( X )= e ~a Xlnin _ e -« max ' Xmin < X < X ma| (25) 

Be-* 

fY(y) = ZgZ > ymin < J/ < Z/max- (26) 

g Pi/mm — g Pymax 

Based on the copula modeling method introduced in Jiang et al. (2009), the coupled probability 
density function fxy{x, y) is 

fxy(x,y) = dF x dF Y 

-0( e -8 - \\ e -8F x {x) -9F Y (y) 

= T ^ ^ ^fx(x)f Y (y), (27) 

e- 9 - 1 + (e-^M*) - l)(e- e ^fe) - 1) 

where the function C(Fx-(x), F Y (y); 9) is given by 



-1 







C(ui, u 2 ; 9) = — In 1 + ^ -f^- i , (28) 



e"»- 1 



Mi is the integral of fx(x), and w 2 is the integral of f Y (y). Thus, 

= F x (x) = — — — — , x mhl < x < x max , (29) 

g " x min — g ux max 
g-pVmin _ e ~/3y 

U 2 = F Y (y) = —r , 2/ m in < V < t/max- (30) 

The dependence parameter (—00 < 9 < 00) can be positive, zero and negative, corresponding 
to the positive dependence, independence and negative dependence between two variables x and 
y, respectively. When 9 approaches to zero, the term d 2 C/dF x dF Y would approach to one, and 
fxy( x , y) = fx{x)f Y (y). Thus, d 2 C/dF x dF Y is called the coupling factor in this paper as it controls 
the x — y dependence. 

The function C(ui,u 2 ;9) in Eq. (28) is called the Frank copula function. In fact, there are 
many available copula functions, and the reason why we choose this one is that it is more flexible 
as it allows to have negative, zero, and positive correlations. 

Further, the expected number of exoplanets in the area dxdy in a survey of N s stars is 

n(x, y)dxdy = cN s f XY (x, y)dxdy, (31) 
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where the parameter c is a constant to be determined. For a given function n(x,y), the number of 
planets with masses between Mi and M 2 and periods between Pi and P2 is determined by 

<-ln(P 2 /P ) r\n(M 2 /M ) 



rln(r 2 /r ) rln(M 2 /M ) 

N[ Mi ,m 2 ][pi,P2] = / n(x,y)dxdy 

J\n(P 1 /P ) Jln(M 1 /M ) 

Hn(P 2 /P ) rln(M 2 /M ) 

= cN s / f XY (x,y)dxdy. (32) 

J\n(Pi/Pn) J\n(Mi/M n ) 



On the other hand, the likelihood function is 



L = Y\n(xi, yi) exp - / n(x,y)dxdy 
i=i L Jd 

= YlcN a fxY(x i ,y i )exp -cN s / / f XY (x,y)dxdy 

, = i L Jymiii Jv+y/3 



(33) 



After the derivation shown in Appendix B, In L is finally expressed as a function of a, j3, and 
9. The maximum likelihood method is used to simultaneously estimate the parameters a, f3 and 9 
through the full log likelihood InL. The estimates of c, a, (3 and 9 for each survey are listed in 
Table 4. Moreover, as the procedure in §4, here we also generalize the result to the case of Multiple 
Surveys and the result is at the bottom of Table 4. The bootstrap method is also used to get error 
bars. Fig. 4 shows the values of a and j3 with error bars. The result of Single Imaginary Survey 
of 175 planets gives a = —0.187 ± 0.034, (3 = —0.133 ± 0.041, so the mass (period) function has 
a power-index —0.813 ± 0.034 (—0.867 ± 0.041). Moreover, the result of Multiple Surveys gives 
a = -0.038 ±0.080, f3 = -0.137 ±0.044, and thus the mass (period) power-index is -0.962 ±0.080 
(—0.863 ± 0.044). On the other hand, the mass (period) power-index obtained in Tabachnik and 
Tremaine (2002) is —1.11 ± 0.10 (—0.73 ± 0.06). Thus, considering the error bars, our results are 
consistent with those in Tabachnik and Tremaine (2002). However, our mean values imply a flatter 
mass function but a slightly steeper period function. We hope to use the future data to investigate 
whether this new result remains to be valid. 

In Fig. 5(a), from the result of Multiple Surveys, the coupled probability density function in 
x — y space, fxv(x,y), is shown as a three-dimensional plot. The corresponding contour is in Fig. 
5(b). To visualize it in a realistic space, the above two are transformed to be in M — P space and 
shown in Figs. 6(a)-(b), where the coupled probability density function in M — P space, /mp(M, P) 
is defined by 

f MP (M,P) = f XY (x,y)\J\, (34) 
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where J is the Jacobian determinant dx/dM x dy/dP. As we know from Eq.(27), the mass- 
period coupling is primarily determined by the coupling factor. In order to visualize it, the three- 
dimensional and color contour plots in x — y space are shown in Figs. 7(a)- (b). It shows that the 
coupling factor is bigger when x and y have the same sign, i.e. both positive or both negative. This 
implies that each group of exoplanets would have its own strength of mass-period correlations if the 
exoplanets are clustered into a few groups. For example, Marchi et al. (2009) did the clustering 
analysis on exoplanets, and found the existence of two types of close-in planets. They discovered 
that these two types of planets have very different distributions of semi-major axes. Moreover, 
Marchi et al. (2009) also proposed possible mechanisms to produce these two different types of 
exoplanets. 

Figs. 8(a)-(b) show the coupled mass-period probability density function (pdf), fxy(x,y), 
for a few given masses or periods. For the purpose of comparison, the pdfs without coupling fac- 
tor, i.e. fx(x) x fy(y) are also plotted in Figs. 8(c)-(d). In Fig. 8(a) and Fig. 8(c), the solid 
curve is for x=\n(lMj/M ), the dotted curve is for x=ln(5Mj/M ), the short dashed curve is for 
x=ln(10Mj/M ), and the long dashed curve is for x=\n(15Mj/M ). In Fig. 8(b) and Fig. 8(d), the 
solid curve is for 2/=ln(lday/P ), the dotted curve is for y=ln(50days/P ), the short dashed curve 
is for y=ln(100days/Po), and the long dashed curve is for y=ln(150days/Po)- in general, it shows 
that the coupling factor does change the shapes of mass and period functions. To be more realistic, 
the curves in Fig. 8 are re-plotted as functions of M and P as shown in Fig. 9. 
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Table 4. 



S 11 FA/PA/ 


Data Set 


N 




8 




r x 10 2 


Lick 


i/v 


7 


f) 531 4- 1 D47 


_n 194 +n S34 


2 289 + Fi 201 


7 cof) + 6 044 


r^nra lip 




38 


n 3i4 + n 158 


-0 S44+ 096 




10 1 01 +1 681 

XW.XWX l 1 .UU1 


r (1 nnip 

LJ1UU1C 




14 

X t: 


_n 4^ + n ^6^ 


-0 203+ 233 


4 7^7 + 2 978 


6 878 + 1 1 234 


HARPS 




23 


046+ 1 Q3 


092+ 1 68 


1 831 + 1 863 

X . OO X _l_ l.OUtJ 


7 01 3 + 312 


N2K 


(E) 


14 


045 ± 438 


130 ± 191 

V 7 > _L '.J V 7 1 V 7 > X r ./ X 


1.795± 2.045 


7 049 ± 1 621 


Keck 


(F) 


48 


0.019±0.132 


-0.217± 0.082 


3.957 ± 1.521 


8.145 ± 0.882 


AAPS 


(G) 


23 


0.018 ± 0.220 


-0.416±0.139 


6.786± 2.288 


7.337 ± 1.353 


Other 


(H) 


17 


-0.069 ± 0.402 


-0.229±0.182 


3.970± 2.511 


7.166 ± 1.635 


Single 




175 


-0.187 ± 0.034 


-0.133± 0.041 


4.848± 0.779 


6.832 ± 0.098 


Multiple 




184 


-0.038 ± 0.080 


-0.137± 0.044 


2.973± 0.639 


7.378 ± 0.289 



7. The Fraction of Stars with Planets 

From both theoretical and observational points of view, it is important to know what fraction 
of stars would have planets. If we divide A^[m 1 ,m 2 ][pi,P2] by N s in Eq. (24), the probability that a 
star could host a planet with masses between Mi and M 2 and periods between Pi and P 2 is given 



a(3 [\M J \M J ] 1 \P J \P 
for the independent mass-period function. 

As the power indexes obtained from the results of Single Imaginary Survey are different from 
those from Multiple Surveys, we calculated both results using Eq. (35). We used the symbol Probs 
to represent the results of Single Imaginary Survey and used Probu for the results of Multiple 
Surveys. 

On the other hand, for coupled mass-period functions, we have 

Probe = — / / n(x,y)dxdy= / / cf X v{x,y)dxdy, 

1S <s J\n(P 1 /P ) Jln(Mi/M ) J J 

-0(p- e - l) e -6Fx(x) -6F Y (y) 

l f x (x)f Y (y)dxdy. (36) 




e- e - 1 + (e-^M*) - lXe-^Ww) - 1)] 



-16- 



Similarly, we used the symbol Probes to represent the results of Single Imaginary Survey and used 
ProbcM for the results of Multiple Surveys in the case of coupled mass-period functions. 

Tabachnik & Tremaine (2002) estimated the expected number of planets per star for given 
period and mass range through their best results of multiple surveys, which is in fact equivalently 
defined as Prob M here. When M 1 = Mj, M 2 = 10Mj, Pi = 2 days, and P 2 = 10 yrs = 3650 days, 
Tabachnik & Tremaine (2002) found Prob M = 0.036 and concluded that 4 percent of solar-type 
stars have a planet in the above ranges. For the same given ranges, we obtain: Probs = 0.02618, 
Prob M = 0.04667, Prob cs = 0.02909, and Prob CM = 0.02273. Thus, our results of Multiple 
Surveys are similar to those in Tabachnik & Tremaine (2002). Moreover, the estimated probabilities 
from the coupled mass-period functions are smaller but still consistent with those in Tabachnik & 
Tremaine (2002). In Tabachnik & Tremaine (2002), the case with M 1 = 0.003M/ (i.e. Earth Mass), 
M 2 = 10Mj, P\ = 2 days, and P 2 = 10 yrs = 3650 days is also estimated and has a probability as 
0.18. However, our results show that: Prob s = 0.06406. Prob M = 0.1264, Probes = 0.06492, and 
ProbcM = 0.0736. Note that the result with Mi = 0.003M/ (i.e. Earth Mass) is an extrapolation 
as this small mass is out of the mass range of 175 samples. Under the assumption that the smaller 
planets shall follow the trend of our mass-period function, and the fact that the Earth-Mass planet 
does exist in our universe, i.e. Solar System, this extrapolation shall lead to a good estimation. 

On the other hand, Naef et al. (2005) also gave estimations about the fractions of star with 
planets more massive than 0.5 Mj within three periods: 0.7 percent for period < 5 days, 4 percent 
for period < 1500 days, and 7.3 percent for period < 3900 days. Correspondingly, using our samples 
and equations here, for Mi = 0.5Mj, M 2 = M max , Pi = P mm , and P 2 = 5 days (Note that, in 175 
samples, the smallest mass M min = 0.0158, the largest mass M max = 18.39, the smallest period 
-Pmin = 1-328, the largest period P max = 5218, and the units are Mj and days), we obtained the 
fraction of stars with planets to be less than 1 percent. For Mi = 0.5Mj, M 2 = M max , Pi = P m i n , 
and P 2 = 1500 days, we obtained the fraction of stars with planets to be about 3 to 6 percent. 
Moreover, for M x = 0.5Mj, M 2 = M max , Pi = P min , and P 2 = 3900 days, the results are about 4 to 
8 percent. 

Finally, we are interested in the possibility to have a planet with mass between Earth Mass 
and Neptune Mass, for any period, so we set Mi = 0.003, M 2 = 0.05, Pi = P m in, and P 2 = P max , 
and get: Prob s = 1.604 %, Prob M = 3.649 %, Prob cs = 1-263 %, and Prob CM = 2.542 %. We are 
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also interested in the possibility between Neptune Mass and Jupiter Mass, for any possible period, 
thus we set Mi = 0.05, M 2 — 1, Pi — P m in, and P 2 = Pmax, and obtain: Probs = 2.589 %, ProbM 
= 5.182 %, Probes = 2.319 %, and Prob C M = 3.021 %. All the above mentioned results calculated 
by the equations in this paper are listed in Table 5. 



Table 5. 



Case 


Mi 


M 2 


Pi 


P 2 


Probs 


Prob M 


Probes 


ProbcM 


1 


1. 


10. 


2. 


3650. 


0.02618 


0.04667 


0.02909 


0.02273 


2 


0.003 


10. 


2. 


3650. 


0.06406 


0.12640 


0.06492 


0.0736 


3 


0.5 


18.39 


1.328 


5. 


0.00453 


0.00790 


0.00154 


0.00169 


4 


0.5 


18.39 


1.328 


1500. 


0.03527 


0.06255 


0.03095 


0.02661 


5 


0.5 


18.39 


1.328 


3900. 


0.04290 


0.07631 


0.04196 


0.03658 


6 


0.003 


0.05 


1.328 


5218.0 


0.01604 


0.03649 


0.01263 


0.02542 


7 


0.05 


1. 


1.328 


5218.0 


0.02589 


0.05182 


0.02319 


0.03021 



8. Concluding Remarks 

In this paper, several steps have been taken to establish the fundamental mass-period functions 
of exoplanets. First of all, by reading a great number of published papers, we constructed a 
reference-based catalog of 175 exoplanets, which were all discovered through the Doppler-shift 
method. Employing this catalog, we determine the independent mass-period functions for individual 
surveys and also for the case of Multiple Surveys. Moreover, the coupled mass-period functions are 
also constructed for both individual and Multiple Surveys. The selection effects of surveys are 
considered in all the results, and thus it is the first time in this field that the selection effect is 
considered in the coupled mass-period functions. Our results are consistent with those in Tabachnik 
and Tremaine (2002) with the main differences that our results imply a flatter mass functions but 
a steeper period function. 

On the other hand, our coupled mass-period functions are used to predict the possible fractions 
of stars in given mass and period ranges. Our results are consistent with previous works, and show 
that about 2.5 percent of stars would have a planet with mass between Earth Mass and Neptune 
Mass, and about 3 percent of stars would have a planet with mass between Neptune Mass and 
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Jupiter Mass. 
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Appendix A 



Data Set 


Survey 


Reference 


Ref. ID 


N 


Planet 


M 


P 


(A) 


Lick 


Marcy et al.(2002) 


(A-l) 


1 


55 Cnc d 


3.835 


5218 






Fischer et al. (2002a) 


(A-2) 


1 


47 Uma c 


1.34 


2594 






Butler et al.(1997) 


(A-3) 


3 


55 Cnc b 


0.824 


14.65162 












Tau Boo b 


3.9 


3.3135 












Ups And b 


0.69 


4.61708 






Johnson et al.(2008) 


(A-4) 


2 


Kappa CrB b 


1.8 


1191 












HD 167042 b 


1.6 


416.1 


(B) 


Coralie 


Mayor et al. (2004) 


(B-l) 


16 


HD 19994 


2 


454 












HD 65216 


1.21 


613.1 












HD 92788 


3.86 


377.7 












HD 111232 


6.8 


1143 












HD 114386 


0.99 


872 












HD 142415 


1.62 


386.3 












HD 147513 


1 


540.4 












HD 196050 


3 


1289 












HD 216437 


2.1 


1294 












HD 216770 


0.65 


118.45 












HD 6434 


0.48 


22.09 












HD 121504 


0.89 


64.6 












HD 83443 b 


0.4 


2.985625 












HD 82943 b 


1.75 


441.2 












HD 82943 c 


2.01 


219 












HD 169830 c 


4.04 


2102 






Tamuz et al. (2008) 


(B-2) 


2 


HD 4113 


1.56 


526.62 












HD 156846 


10.45 


359.51 






Udry et al.(2003) 


(B-3) 


1 


HD 73256 


1.87 


2.54858 






Naef et al.(2001a) 


(B-4) 


3 


GJ 3021 


3.32 


133.82 












HD 52265 


1.13 


118.96 












HD 169830 b 


2.88 


225.62 






Queloz et al.(2000) 


(B-5) 


1 


GJ 86 


4.01 


15.766 
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burvey 


Keterence 


Ket. ID 


AT 
IN 


Planet 


71 AT 

M 


J ) 

r 


/t->\ 

(B) 


Corahe 


tt 1 x i/nnnn\ 

Udry et al(2000) 


(B-6) 


2 


t t t ~\ ^rroon 

HD 75289 


0.42 


3.51 












HD 130322 


A f\ o 

1.08 


10.724 






Santos et al. (21)00) 


(B-7) 


1 


HD 192263 


0.72 


24.348 






Ox x 1 /nnm \ 

Santos et al.(200l) 


(B-8) 


2 


HD 28185 


5.7 


o o o 

383 












tttv m on /in 

HD 213240 


4.5 


951 






Correia et al.(2005) 


(B-9) 


1 


HD 202206 c 


2.44 


1383.4 






Pepe et al. (2002) 


(B-10) 


2 


HD 108147 


n ,1 

0.4 


1 n nn -i 

10.901 












HD 168746 


n o O 

0.23 


6.403 






TT1 x 1 / C\C\C\C\\ 

Udry et al.(2002) 


(B-ll) 


4 


T T T ~\ AAA /"I O ^7 

HD 141937 


n T 

9.7 


653.22 












HD 162020 


13.75 


o ,moi no 

8.428198 












HD 168443 c 


1 o 1 

18.1 


1765.8 












HD 202206 b 


17.4 


255.87 






Zucker et al. (2004) 


/ID 1 0\ 

(B-12) 


1 


HD 41004 A 


2.3 


655 






oantos et al. (2002) 


(B-13) 


1 


HD 41UU4 B 


in /) 


1.3283 






bggenberger et al.(200b) 


(B-14) 


1 


HD 14zUzz 


4.4 


1923 


(C) 


Elodie 


r \ 1 1 lxl /r>nnr\ 

Galland et al. (2005) 


(C-l) 


1 


HD 33564 


n 1 

9.1 


o O O 

388 






Perrier et al. (2003) 


(C-2) 


6 


T T T ~\ o f ry a 

HD 8574 


2.23 


228.8 












HD 23596 


7.19 


1558 












HD 33636 


n o o 

9.28 


2127.7 












HD 50554 


/i n 

4.9 


1 oTn 

1279 












HD 106252 


6. 81 


1 r - nn 

1500 












HD 190228 


/i nn 

4.99 


1127 






at £ x l /rinn >i \ 

Naei et al.(2004) 


(C-3) 


3 


HD 74156 b 


1.88 


51.65 












T T T ~\ ry a 1 r - z" 1 

HD 74156 c 


8.03 


2476 












t t t ~\ 1 /i rr"7r / 1 /i tt \ 

HD 145675(14 Her) 


4.64 


1773.4 






iNaet et al.(2U0o ) 


(C-4) 


1 


tiD lyUooO b 


1.002 


ocn 1 






Moutou et al.(2006) 


(C-5) 


1 


HD 185269 


0.94 


6.838 






da. Silva et al.(2006) 


(C-6) 


1 


HD 118203 


2.13 


6.1335 






Mayor & Queloz (1995) 


(C-7) 


1 


51 Peg 


0.468 


4.23077 


(D) 


HARPS 


Udry et al. (2006) 


(D-l) 


1 


HD 4308 


0.047 


15.56 
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Data Set 


Survey 


Reference 


Ref. ID 


N 


Planet 


M 


P 


(D) 


HARPS 


Moutou et al.(2005) 


(D-2) 


3 


HD 2638 


0.48 


3.4442 












HD 27894 


0.62 


17.991 












HD 63454 


0.38 


2.81782 






Pepe et al. (2004) 


(D-3) 


1 


HD 330075 b 


0.76 


3.369 






Lo Curto et al.(2006) 


(D-4) 


1 


HD 212301 


0.45 


2.457 






Lovis et al.(2005) 


(D-5) 


3 


HD 93083 


0.37 


143.58 












HD 10193 


0.3 


70.46 












HD 102117 


0.172 


20.67 






Santos et al.(2004) 


(D-6) 


f 


HD 160691 d 


0.044 


9.55 






Pepe et al.(2007) 


(D-7) 


f 


HD 160691 e 


0.5219 


310.55 






Bonfils et al.(2005) 


(D-8) 


1 


Gl 581 b 


0.0492 


5.3683 






Udry et al.(2007) 


(D-9) 


2 


Gl 581 c 


0.0158 


12.932 












Gl 581 d 


0.0243 


83.6 






Bonfils et al.(2007) 


(D-fO) 


f 


GJ 674 


0.037 


4.6938 






Melo et al.(2007) 


(D-ll) 


f 


HD 219828 


0.066 


3.8335 






Santos et al. (2007) 


(D-12) 


f 


HD 171028 


1.83 


538 






Naef et al. (2007) 


(D-13) 


3 


HD 100777 


1.16 


383.7 












HD 190647 


1.9 


1038.1 












HD 221287 


3.09 


456.1 






Lovis et al. (2006) 


(D-14) 


3 


HD 69830 b 


0.033 


8.667 












HD 69830 c 


0.038 


31.56 












HD 69830 d 


0.058 


197 


(E) 


N2K 


Wright et al. (2007) 


(E-l) 


3 


HIP 14810 b 


3.84 


6.6742 












HIP 14810 c 


0.76 


95.2914 












HD 154345 b 


2.03 


10900 






Johnson et al. (2006) 


(E-2) 


3 


HD 33283 


0.33 


18.179 












HD 86081 


1.5 


2.1375 














n 7i 


Z,U . t o 






Fischer et al.(2006) 


(E-3) 


2 


HD 149143 


1.33 


4.072 












HD 109749 


0.28 


5.24 






Fischer et al.(2007) 


(EM) 


5 


HD 11506 


4.85 


1280 














Q O 
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Data Set 


Survey 


Reference 


Ref. ID 


N 


Planet 


M 


P 


(E) 


N2K 








HD 170469 


0.67 


1145 












HD 17156 b 


3.111 


21.21725 






Sato et al.(2005) 


(E-5) 


1 


HD 149026 


0.36 


2.8766 






Fischer et al.(2005) 


(E-6) 


1 


HD 88133 


0.22 


3.41 


(F) 


Keck 


Butler et al. (2003) 


(F-l) 


6 


HD 108874 b 


1.36 


395.4 












HD 114729 


0.82 


1131.478 












HD 72659 


2.96 


3177.4 












HD 128311 b 


2.18 


448.6 












HD 49674 


0.115 


4.9437 












HD 37124 c 


0.683 


2295 






Marcy et al.(2005) 


(F-2) 


5 


HD 183263 


3.69 


634.23 












HD 117207 


2.06 


2627.08 












HD 188015 


1.26 


456.46 












HD 45350 


1.79 


890.76 












HD 99492 


0.109 


17.0431 






Robinson et al.(2007) 


(F-3) 


2 


HD 75898 b 


1.48 


204.2 












HD 5319 b 


1.94 


675 






Butler et al.(1998) 


(F-4) 


1 


HD 187123 b 


0.52 


3.097 






Marcy et al.(1999) 


(F-5) 


2 


HD 210277 


1.23 


442.1 












HD 168443 b 


8.02 


58.11289 






Johnson et al.(2007) 


(F-6) 


1 


GJ 317 


1.2 


692.9 






Vogt et al.(2000) 


(F-7) 


6 


HD 10697 


6.12 


1077.906 












HD 37124 b 


0.61 


154.46 












HD 134987 


1.58 


260 












HD 177830 


1.28 


391 












HD 192263 


0.72 


24.348 












HD 22258 


5.11 


572 










i 

L 


vjr J O^j 




lOc/U 






Vogt et al. (2002) 


(F-9) 


5 


HD 4203 


1.65 


400.944 












HD 4208 


0.8 


812.197 












HD 33636 


9.28 


2127.7 














1 n 


a o7<; 
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Data Set 


Survey 


Reference 


Ref. ID 


N 


Planet 


M 


P 


(F) 


Keck 


Marcy et al.(2001a) 


(F-10) 


1 


HD 168443 c 


18.1 


1765.8 






Rivera et al. (2005) 


(F-ll) 


1 


GJ 876 d 


0.018 


1.93776 






Vogt et al.(2005) 


(F-12) 


7 


HD 128311 c 


3.21 


919 












HD 50449 


1.71 


2582.7 












HD 37124 d 


0.6 


843.6 












HD 190360 c 


0.057 


17.1 












HD 108874c 


1.018 


1605.8 












HD 37124 c 


0.683 


2295 












HD 217107 c 


2.5 


3352 






Butler et al.(2004) 


(F-13) 


1 


GJ 436 


0.072 


2.64385 






Butler et al. (2006a) 


(F-14) 


5 


HD 11964 b 


0.11 


37.82 












HD 66428 


2.82 


1973 












HD 99109 


0.502 


439.3 












HD 107148 


0.21 


48.056 












HD 164922 


0.36 


1155 






Marcy et al.(2000) 


(F-15) 


2 


HD 16141 


0.23 


75.56 












HD 46375 


0.249 


3.024 






Butler et al.(2000) 


(F-16) 


2 


HD 52265 


1.13 


118.96 












BD 103166 


0.48 


3.488 


(G) 


AAPS 


Tinney et al.(2003) 


(G-l) 


4 


HD 73526 b 


2.9 


188.3 












HD 76700 


0.197 


3.971 












HD 30177 


9.17 


2819.654 












HD 2039 


4.85 


1192.582 






Jones et al. (2002b) 


(G-2) 


3 


HD 196050 


3 


1289 












HD 216437 


2.1 


1294 












HD 160691 c 


3.1 


2986 






Jones et al. (2002a) 


(G-3) 


1 


HD 39091 


10.35 


2063.818 






Tinney et al.(2002) 


(G-4) 


2 


HD 142 


1 


337.112 












HD 23079 


2.61 


738.459 






Jones et al.(2006) 


(G-5) 


2 


HD 187085 


0.75 


986 












HD 20782 


1.8 


585.86 



- 29 - 



Data Set 


Survey 


Reference 


Ref. ID 


N 


Planet 


M 


P 


(G) 


AAPS 


O'Toole et al.(2007) 


(G-6) 


2 


HD 23127 


1.5 


1214 












HD 159868 


1.7 


986 






Butler et al.(2001) 


(G-7) 


2 


HD 160691 b 


1.67 


654.5 












HD 27442 


1.28 


423.841 






Carter et al.(2003) 


(G-8) 


1 


HD 70642 


2 


2231 






Tinney et al.(2001) 


(G-9) 


1 


HD 179949 


0.95 


3.0925 






Jones et al.(2003) 


(G-10) 


1 


Tau Gruis b 


1.49 


1442.919 






Tinney et al.(2006) 


(G-ll) 


1 


HD73526 c 


1.6 


416.1 






Tinney et al.(2005) 


(G-12) 


3 


HD 117618 


0.19 


52.2 












HD 208487 


0.45 


123 












HD 102117 


0.172 


20.67 


(H) 


Others 


Fischer et al. (2001) 


(H-l) 


3 


HD 12661 b 


2.3 


263.6 












HD 92788 


3.86 


377.7 












HD 38529 b 


0.78 


14.30 






Marcy et al. (2001b) 


(H-2) 


1 


GJ 876 c 


0.56 


30.1 






Fischer et al. (2003) 


(H-3) 


3 


HD 40979 


3.32 


267.2 












HD 12661 c 


1.57 


1444.5 












HD 38529 c 


12.7 


2174.3 






Delfosse et al.(1998) 


(H-4) 


1 


GJ 876 b 


1.935 


60.94 






Cochran et al. (1997) 


(H-5) 


1 


16 Cyg B b 


1.68 


799.5 






Sozzetti et al. (2006) 


(H-6) 


1 


HD 81040 


6.86 


1001.7 






Naef et al. (2001b) 


(H-7) 


1 


HD 80606 b 


3.41 


111.78 






Fischer et al. (1999) 


(H-8) 


2 


HD 195019 A b 


3.7 


18.20163 












HD 217107 


1.33 


7.12689 






Butler et al.(1999) 


(H-9) 


2 


Ups And c 


1.98 


241.52 












Ups And d et al. 


3.95 


1274.6 






Fischer et al. (2002b) 


(H-10) 


1 


HD 136118 


11.9 


1209 






Korzennik et al.(2000) 


(H-ll) 


1 


HD 89744 


7.99 


256.605 
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Appendix B 



The derivation of the log likelihood InL of coupled mass-period functions is shown here. As in 
Tabachnik and Tremaine (2002), the likelihood function is 



A' 



L = Y\n(xi,yi) exp - / n(x,y)dxdy 
i=i L Jr 

N 

= \\cN 8 fxY(xi,yi) exp 



»=i 



cN x I / f XY (x,y)dxdy 

2/min Jv+y/3 



(37) 



By the other equations in §6, we have 

r [ x — t r f Xn ^ d 2 C(F x {x),F Y {y);6) 



/ fxY(x,y)dxdy -- 

'Vmin Jv+y/3 Jy min Jv+y/3 Ot X Ot Y 



fx(x)f Y (y)dxdy, 



, fi | ^(Fx(^max),Fy(l/)^) _ aC(Fy(t; + y/3),Fy(y);g) - . , ^ . 

ymin V 



•{ I 9Fy cVfy 

= C(F x (:r max ), Fy(u); 0) - C(F x (x max ), Fy(y min ); 0) 



dC(F x (v + y/3\F Y (y);9) 



2/min 



r e -0(Fy(j/)+F x (^+s,/3)) _ e -0Fy(j/)l 



e -e + e -0(Fy( J /)+F x (^+ J //3)) _ e -eF Y {y) _ e -eF x (v+y/3) 



^fr(y)dy, 



fv(y)dy, 
(38) 



where C(F x (x max ), Fy(</ min ); 9) = C(1,O,0) = ,and C(F x (x max ), 0) = C(l,F y («);0) = 

F y (w). From Eq.(37) and Eq.(38), we have 



n(x, y)dxdy 



D 



exp 

= 6XPl ' " " '• ' 1 " ' ' e -0 _|_ e -e(Fy(j/)+Fx(«+y/3)) _ e -9F Y (y) _ e -9F x (v+y/3) ' 

Vmin 



\-cn s (f Y (u)- r 

< J 2/mii 



Fy(7, 

and the log-likelihood is 



[ e -9(F y (v)+F x (t;+i//3)) _ e" eFy(?/) ]/y (y) 



dy) L (39) 



N N 

InL = N]n(cN s ) - a^ Xi - fi^yi + iVln (- 

i=i i=i ' 



Ci;f' m j n p ftx max 



) +iVln (- 







-/3j/, 



mm — p max 



N 

+ J] In (C ttlU2 (F Y (x i ), F Y {y t ); 0)) - ciV s /(a, /3, u, v), 



(40) 



where 



C UlU2 (F x (x i ),F Y (y i );< 



-Q(er e - i) e ~ eF x{^) e -eF Y { yi ) 



e- e - 1 + (e-^(^) - l)(e-^(w) - 1)] 



2 ' 



I(a, (3,u,v) = F Y (u 



J Vmu 



[ e -9(F Y (y)+F x (v+y/3)) _ e - 9FY ^]f Y (y) 
e -9 _|_ e -e(Fy(»)+F x (t)+j//3)) _ e -6F Y {y) _ e -eF x (v+y/3)) 



(41) 



dy. (42) 
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By dlnL/dc = 0, we obtain c = N/[N s I(a, (3, u, v )]. It is then substituted into Eq.(40) and we 

have 

N N 

InL = N In(iV) - N ln(J) - a ^ x % - (5 £ + N In ( g ) 

i=l i=l 

+An » _ e -/w ) + E ln *Hk); *)) ~ N. (43) 

1=1 
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Fig. 1. — The estimators of the exponents, a and (3, with error bars for all surveys listed in Table 
1, and also the results of Multiple Surveys. 
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Fig. 2. — The distributions of exoplanets in M and P space. The solid lines are the histograms 
of 175 exoplanets. By Eq.(24), the crosses with dotted curves are the results of Single Imaginary 
Survey, and the circles with dashed curves are the results of Multiple Surveys. 
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Fig. 4. — The estimators of the exponents, a and (3, with error bars for all surveys listed in Table 
4. 
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r. -c 




Fig. 5. — (a) The coupled probability density function in x — y space, i.e. fxy{x, y), from the results 
of Multiple Surveys, (b) The corresponding contour plot. 
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0.5 1 1.5 2 2.5 5 3.5 4 4.5 5 5.5 

M 



Fig. 6. — (a) The coupled probability density function m M — P space, i.e. fMp{M,P), from the 
results of Multiple Surveys, (b) The corresponding contour plot. 




Fig. 7. — (a) The three-dimensional plot of the coupling factor for the results of Multiple Surveys, 
(b) The corresponding contour plot. 
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Fig. 8. — The mass and period functions in x — y space for the results of Multiple Sur- 
veys, (a) The period functions of x=ln(lMj/M ) (solid curve), x=\n(5Mj/M ) (dotted curve), 
x=ln(10Mj/M ) (short dashed curve), and x=ln(15Mj/M ) (long dashed curve), (b) The mass 
functions of |/=ln(lday/Po) (solid curve), y=ln(50days/Po) (dotted curve), t/=ln(100days/Po) (short 
dashed curve), and ?/=ln(150days/P ) (long dashed curve), (c) The independent period func- 
tions of x=\n(lMj/M ) (solid curve), x=\n(5Mj/M ) (dotted curve), x=ln(10M//M ) (short 
dashed curve), and x=ln(15Mj/M ) (long dashed curve), (d) The independent mass functions 
of y=ln(lday/P ) (solid curve), y=ln(50days/P ) (dotted curve), |/=ln(100days/P ) (short dashed 
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Fig. 9. — The mass and period functions in M — P space for the results of Multiple Surveys, (a) The 
period functions of M — lMj (solid curve), M = 5Mj (dotted curve), M = 10Mj (short dashed 
curve), and M = 15Mj (long dashed curve), (b) The mass functions of P = 1 day (solid curve), 
P = 50 days (dotted curve), P = 100 days (short dashed curve), and P = 150 days (long dashed 
curve), (c) The independent period functions of M — lMj (solid curve), M = 5Mj (dotted curve), 
M = 10Mj (short dashed curve), and M = 15Mj (long dashed curve), (d) The independent mass 
functions of P = 1 day (solid curve), P = 50 days (dotted curve), P = 100 days (short dashed 
curve), and P = 150 days (long dashed curve). 



